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Abstract 

The changes of polarizational, elastic and structural properties of mutiferroic BiFeC>3 under [111] 
direction uniaxial stress are calculated using density functional theory within the Perdew-Burke- 
Ernzerhof revised for solids (PBEsol) + U approach, and compared with available measurements 
or predictions. A large ferroelectric polarization of 87.5 /iC/cm 2 is found, agreeing with other 
theoretical and experimental values, and the polarization changes near-linearly within the uniaxial 
stress range of 8 GPa to -8 GPa. This property can be used to indirectly control the radiative 
recombination of luminous organ and the electrical properties of transistor structure. In addition, 
we have also investigated the elastic stiffness coefficients and the elastic compliance coefficients in 
the same uniaxial stress range, which provide helpful guidance for the applications of BiFeC>3. 

PACS numbers: 77.84.-s, 77.22.Ej, 62.20.D- 
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I. INTRODUCTION 



Multiferroics are promising materials because of their fundamental physical properties as 
well as their wide range of potential applications^-. Among these materials, BiFeOa is the 
most promising multiferroics from the applications perspective^, and the latest results, such 
as photovoltaic effect^ and ferroelectric-resistive random access memory^, it offers continue 
to fuel the interest in it. 

At room temperature, bulk BiFeC>3 has long been known to be ferroelectric, with a Curie 
temperature ~1110 K, and antiferromagnetic, with a Neel temperature ~643 K 6 ~ 8 . It is a 
rhombohedral R3c perovskite with a large electric polarization 9 88.7 fiC/cm 2 pointing along 
the [111] direction, and Bi, Fe, and O are displaced relative to one another along this 
threefold axis. The largest relative displacements are those of Bi relative to O, consistent 
with a stereochemical^ active Bi lone pair which induces a structural distortion that can 
lead to ferroelectricity 9 . But under high pressure (i.e. strain), the Bi atom tend to form 
planar BiC>3 groups and the Fe 3+ show a high spin (HS) - low spin (LS) crossover transition 
at pressures around 40 GPa^. Furthermore, when the rhombohedral perovskite BiFeC>3 
is grown on a substrate having a square in-plane lattice, its symmetry is lowered to a 
monoclinic phase 11 . As the epitaxial strain increaseing, it undergoes a phase sequence of 
rhombohedral (R) to monoclinic (i?-like Ma) to monoclinic (T-like Mq) to tetragonal (T), 
which is otherwise seen only near morphotropic phase boundaries (MPB) in lead-based 
solid-solution perovskites 12 . The above results thus show that "strain engineering" is an 
important tool both in fundamental studies to understand complex phase transition as well 
as in approaches to find new, lead-free materials with technologically relevant properties, 
such as large piezoelectric coefficients and tunable polarization. In this work, we report the 
changes of polarization and elasticity properties under uniaxial stress (strain) in BiFe03. 

Elastic strain is an effective way to control the magnetic and the polarization of 
BiFeO^-^r— . Ramazanoglu et air showed that small uniaxial pressure (i.e. elastic strain) 
will produce significant changes in the populations of magnetic domains and rotate the mag- 
netic easy plane of domains in BiFeCV Jang et al. 13 provided direct experimental evidence 
that (OOl)p-oriented BiFeC>3 films exhibit a strong strain tunability of their out-of-plane 
remanent polarization. However, no theoretical investigation has been done until now to 
study the changes of polarization and elasticity properties under uniaxial stress in BiFeC>3, 
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which inhibits the applications of BiFeOs not only in bulk materials but also in thin films. 
This is another motivation for pursuing this work. 

In this work we analyze the changes of polarization, elasticity and structure properties 
of BiFeC>3 under [111] direction uniaxial stress and discuss its physical mechanism. It is 
of great interest to note that the polarization of BiFeOs changes near-linearly in the stress 
range of 8 GPa to -8 GPa. This property is conducive to the application of BiFeOs, such as 
by controlling the spontaneous polarization to indirectly control the radiant recombination 
of luminous organ and the electrical properties of transistor structurei^^. Furthermore, 
the investigation of elastic and structural properties reveals that the structure of BiFeOs 
is stable within the uniaxial stress range of 8 GPa to -8 GPa, which provides direction for 
future applications. 

II. METHODOLOGY 

At room temperature, BiFeOs is a rhombohedral (space group R3c) G-type antiferro- 
magnetic (AFM) perovskite structure^^. In order to exert stress on [111] direction expe- 
diently, we do not study the rhombohedral cell (a = (3 = 7 7^ 90°) directly but its equivalent 
one, i.e., the hexagonal cell (a = (3 = 90°, 7 = 120°) (see Fig. [1]^. The [001] direction of 
hexagonal cell is the [111] direction of rhombohedral cell. For the sake of calculating the 
uniaxial stress along the [001] direction, we apply a strain in the [001] direction, and then 
keep the lattice constant of [001] direction unchanged and relax its perpendicular lattice 
vectors and all the internal atomic positions. The optimization is done until the other two 
components of the stress tensor (i.e. an and 022) are all smaller than 0.05 GPa. As a result, 
with the strain changing continuously we get the variational uniaxial stress. 

To calculate the structure, polarization, and elastic coefficient of BiFeOs, we use density 
functional theory within the generalized gradient approximation (GGA)^ and GGA+£7 
method^ (with the Perdew-Burke-Ernzerhof revised for solids (PBEsol)^ 1 ^) as implemented 
in the VASP package^. The PBEsol+£7 method was developed by Dudarev et al.—, which 
depended on the effective Coulomb interaction U e ff = U — J. In the present work, U e ff = 4 
eV is used, as it leads to qualitatively and semiquantitatively correct results for all the 
properties investigated^ 1 ^. All results are obtained using the projector augmented wave 
method to represent the ionic cores 2 ^, solving for the following electrons: Fe's 3p 6 3d 6 4s 2 ; 
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Bi's 5d 10 6s 2 6p 3 ; and O's 2s 2 2p 4 . In the first-principles calculation, we employ a plane- 
wave cutoff of 500 eV and k-point sampling of 4x4x2 (comparing with k-point sampling 
of 4x4x8, the error of polarization is less than 0.1%). The electronic energy is convergent 
to at least 0.003 meV/atom, and all force components are relaxed to at least 2 me v/A. 
The optimizations of cell shape and ion's position are performed by Gaussian smearing 
technique 2 ^. The polarization and the elastic coefficients are calculated by Berry phase 2 ^ 1 ^ 
and the finite difference method 2 -^ respectively. 

III. RESULTS AND DISCUSSIONS 

A. Properties of unstrained BiFeC-3 

First of all, we study the structure, polarization and the elastic stiffness coefficients of 
unstrained BiFeO"3, which are all listed in Table HI An earlier report^ reveals that the ferro- 
electric properties of BiFeC>3 are very sensitive to the small changes in the lattice parameters. 
Consequently, for predicting the ferroelectric properties correctly it is needful to estimate 
the structural parameters of BiFeC>3 accurately. As can be seen in Table [J our optimized 
lattice parameters are in good agreement with the experimental values^ and other theoret- 
ical values^, which guarantees the veracity of predicted ferroelectric properties and elastic 
properties. The polarization (87.5 /iC/cm 2 ) is a little bit smaller than the experimental 
values&2£ but almost equal to the theoretical values (88.7 fiC/cm 2 ) in Ref. f2]. The elastic 
stiffness coefficients, except for C44, are all slightly larger than the values in Ref. Taking 
one with another, all the above properties of unstrained BiFeC>3 are consistent with recent 
studies**^. 

B. Near-linearly variational polarization 

The changes of polarization during the stress changing from 16 GPa (tensile stress) to 
-52 GPa (compressive stress) are showed in Fig. [2J Interestingly we observed that the 
polarization P increases almost linearly from 75.45 /zC/cm 2 to 97.08 /iC/cm 2 during the stress 
changing from 8 GPa to -8 GPa. The changing trend agrees with the results in Ref. [2i], 
which investigated the change of polarization in ±3% strained BiFe03, but is opposite to the 
results in Ref. p^], which studied the change of polarization under uniaxial stress in PbTi03. 
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An earlier report^ reveals that the polarization contribution arising from displacement of Fe 
and O atoms are almost canceling out and more than 98% of the net polarization present in 
BiFeOs is contributed by the Bi ions. Another report^ indicates that the polarization comes 
from the lone electron pair of Bi 3+ ions. So, it is reasonable to deduce that the near-linearly 
variational polarization is related to (i) the displacement and (ii) the Born effective charge of 
Bi ions. And it is indeed interesting to analyse the relative displacement and Born effective 
charge of Bi ions in BiFeOs i n detail under diferent uniaxial stress. 

Let us now concentrate on the relative displacement (Am = u - uq) of ions in the cell 
(under fractional coordinates) vs. uniaxial stress. Strikingly, Fig. E^a) reveals that the 
changing curve of relative displacement of Bi 3+ ions is almost perfect linear in the uniaxial 
stress range of 8 GPa to -8 GPa and the moving direction of Bi 3+ ions happen to be [001] 
direction (polar axis, not shown), which are compatible with the near-linearly variational 
polarization. Note that the curves of Fe 3+ ion and 2 ~ ion are almost overlap in the uniaxial 
stress range of 8 GPa to -8 GPa ,which indicates that the Fe06 octahedron is not destroyed 
in this stress range. In addition, not only Bi 3+ ions, but also Fe 3+ and O 2- ions have 
near-linearly changing curve of relative displacement in the uniaxial stress range of 8 GPa 
to -8 GPa. But the curve-slopes of Fe 3+ and O 2- ions are far less than Bi 3+ ions', and 
the changing direction of relative displacement between Bi 3+ ions and Fe 3+ (or O 2- ) ions 
in the cell is opposite (because of opposite slopes), which implys that the largest relative 
displacements are those of Bi 3+ ions relative to 2 ~ (or Fe 3+ ) ions, as revealed in most 
reports of unstrained BiFeOs 3 -''^. 

To get more detailed information on near-linearly variational polarization, the Born ef- 
fective charge (Z*) of ions in BiFe03 were obtained in the stress range of 16 GPa to -52 GPa 
, shown in Fig. [3](b). Interestingly, the |Z*| (absolute value of Z*) of Fe 3+ and 2 ~ ions 
decrease synchronously with the uniaxial stress changing from 8 GPa to -8 GPa, associating 
with the results of relative displacement of Fe 3+ and O 2- ions demonstrated in Fig. [3](a), 
which show that the polarization contribution coming from Fe and O atoms are still almost 
canceling out and the net polarization present in BiFeOs is mainly contributed by the Bi 
ions. However, it should be noted that the Z* of Bi 3+ ions decrease first and then increase, 
not always increase, with the uniaxial stress changing from 8 GPa to -8 GPa. So, although 
the Bi 3+ ions have a biggish relative displacement in the range of 8 GPa to -8 GPa, the 
polarization of BiFeOs increases slowly. This result is in line with the earlier reports 3 -^ 2 - 
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and give a possible explanation for it. 

C. Elastic properties under uniaxial stress 

The elastic stiffness coefficients (c^'s) and elastic compliance coefficients (sy's) as a func- 
tion of uniaxial stress (-033) are shown in Figs. H^a) and|4^b). The cy's and s^-'s collectively 
referred to as elasticity coefficient. They represent whether a material is easy to be stretched 
or not. The bigger the c^-'s, the stronger the resisting deformation ability. On the contrary, 
the bigger the s^'s, the easier a material to be stretched. According to the symmetry of 
R3c space group, the c^'s only have six independent variables 33 , namely, en, en , C13, Cu, 
C33, C44 (The Si/s have six similar variables, viz. sn, S12 , S13, su, S33, S44). All of them can 
be obtained through the relationship between strain and stress. It must be noted that both 
Cjj's and s^-'s should satisfy the Born stability criterions. If not, the structure is unstable 3 - 3 -. 
For a rhombohedral structure, the Born stability criterions are 

cu - |ci 2 | > 0, (I) 

(en + c 12 )c 33 - 2cj 3 > 0, (2) 

(cu - ci 2 )c 44 - 2c 2 14 > 0. (3) 

The form of constraint condition of sy's is the same as the above formulas (not shown). 
Figure BJa) shows that, except for C14, the other e^'s increase near-linearly when the stress 
change from 8 GPa to -8 GPa. When the tensile stress is bigger than ff GPa, there is 
a great change in c^-'s: except for c u , the other c^-'s decrease rapidly and tend to 0. It 
is obvious that the structure is mechanically unstable at larger tensile stress based on the 
Born stability criterions Eq. ([IJ-Q. For example, according to our calculation, the elastic 
stiffness coefficients are cu = 185 GPa, c\i = 104 GPa, C13 = 33 GPa, C14 = 22 GPa, C33 
= 135 GPa and C44 = 11 GPa when the tensile stress is 11.12 GPa. These coefficients can 
not meet the Eq. So the structure is predicted to be elastic instability, due to the 

decrease in c 44 and the increase in c u . This predicted instability could suggest a ferroelastic 
phase transition 3 - or melting 3 ^ achieved at negative pressures (such as in epitaxial films) 
or high temperatures. These results are consistent with an earlier report of BiFeOs under 
hydrostatic pressure^. The difference is that all the q/s of AFM BiFeOs are very sensitive 
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to the hydrostatic pressure (see Fig . 2 of Ref.-). But when we focused on the c^'s in 
Fig. Eta), we only found C33 increases unstably and relative-rapidly due to the high uniaxial 
compressive stress in the [001] direction. 

Similarly, Fig. HJb) shows that, except for S44, the absolute values of s^-'s decrease near- 
linearly during the stress changing from 8 GPa to -8 GPa. When the tensile stress is bigger 
than 11 GPa, the absolute values of su, S12, su and S44 increase rapidly. These great change 
make s^-'s no longer satisfy the Born stability criterions^ 3 - and imply that the structure is 
elastic instability, which is tally with the above analysis. 

D. Structural properties under uniaxial stress 

From the analysis of elastic properties of BiFeOs under uniaxial stress, we found that the 
structure of BiFeOs is stable within the stress range of 8 GPa to -8 GPa, and when the tensile 
stress is greater than 11 GPa, the structure of BiFeOs becomes unstable obviously. This 
predicted structural stability from the analysis of elastic properties should be reflected in 
the change of structural parameters. To validate this idea, we analyzed the relative volume 
V/Vo and the strains e\ = (a/ao — 1), £2 = (c/co — 1) as a function of uniaxial stress, shown 
in Fig. [5j It is obvious that the V/Vq, E\ and e<i change near-linearly during the stress 
changing from 8 GPa to -8 GPa. Their slopes are -4.64 x 10" 3 , 1.67 x 10~ 3 and -7.16 x 10~ 3 , 
respectively. These phenomena indicate that the structure is stable in this uniaxial stress 
range. However, the relationship between V/Vq (or £2) and uniaxial stress becomes nonlinear 
as the tensile stress increases, implying that the strain has crossed the elastic limit (or has a 
ferroelastic phase transition 34 ) and the structure is mechanically unstable. The results are 
consistent with the above analysis of elasticity and provide direction for future applications. 

IV. CONCLUSIONS 

We have studied the changes of polarization, elastic properties and structural stabil- 
ity of BiFeOs with rhombohedral G-type AFM structure under uniaxial stress. First, we 
compare our results with available measurements and predictions. Then, we observe that 
the polarization P increases almost linearly from 75.45 /iC/cm 2 to 97.08 yuC/cm 2 during the 
stress changing from 8 GPa to -8 GPa and we have analyzed its two possible causes: (i) 



7 



the displacement and (ii) the Born effective charge of Bi ions. Third, we study the elastic 
properties and structural stability comprehensively and find that the structure is stable in 
the stress range of 8 GPa to -8 GPa. These investigates provide helpful guidance for the 
applications of BiFeC>3, and will lead to more extensive studies, such as the changing rules 
of elastic properties, magnetic properties, polarization and phase transition under biaxial 
stress in BiFeOs thin films, and so on. 
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TABLE I: The lattice parameters, polarization and the elastic stiffness coefficients of unstrained 
BiFeOs: our calculated results, theoretical value in Ref. f^^jand experimental values. The units 
of a (or c), P, and cy's are angstrom (A), ^C/cm 2 and GPa, respectively. 
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c 44 
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164 


41 


Theory 


5.51 a 


13.54 a 


88.7793.3° 


222 c 


110 c 


50 c 


16 c 


150 c 


49 c 


Experiment 


5.57 d 


13.86 d 


59-90 e /100^ 















Tirst-principles prediction by LSDA+C/, U eff = 4 eV (Refi22). 

h First- principles prediction by GGA (Refj^). 

c First- principles prediction by GGA+C, Ueff = 6 eV (Ref. 7 ). 

d Reference2&. 

e Thin films (400 - 100 nm) on SrRu0 3 /SrTi0 3 (Ref 4). 
-'Reference^. 




FIG. 1: (Color online) Structure of BiFeO"3 with hexagonal cell. 
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FIG. 2: (Color online) The polarization as a function of uniaxial stress -033. 
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FIG. 3: (Color online) (a) Relative displacement (Am = u — uq) and (b) Born effective charge of 
ions in the cell as a function of uniaxial stress -033, where uq is the position of ions in unstrained 
BiFeOa under fractional coordinates. 
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FIG. 4: (Color online) (a) Elastic stiffness coefficients cw's and (b) elastic compliance coefficients 
Sij 's as a function of uniaxial stress -033 . 
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FIG. 5: (Color online) (a) V/Vo, (b) £\ = (a/ao — 1) and £2 = (c/co — 1) as a function of uniaxial 
stress -CJ33, where Vq and ao, cq are the volume and lattice parameters of unstrained BiFeOs. 
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